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ABSTRACT 

Large-scale asymmetries in the stellar mass distribution in galaxies are be- 
lieved to trace non-equilibrium situations in the luminous and/or dark matter 
component. These may arise in the aftermath of events like mergers, accretion, 
and tidal interactions. These events are key in the evolution of galaxies. In 
this paper we quantify the large-scale lopsidedness of light distributions in 25155 
galaxies at 2; < 0.06 from the Sloan Digital Sky Survey Data Release 4 using the 
171 = 1 azimuthal Fourier mode. We show that the lopsided distribution of light 
is primarily due to a corresponding lopsidedness in the stellar mass distribution. 
Observational effects, such as seeing, Poisson noise, and inclination, introduce 
only small errors in lopsidedness for the majority of this sample. We find that 
lopsidedness correlates strongly with other basic galaxy structural parameters: 
galaxies with low concentration, stellar mass, and stellar surface mass density 
tend to be lopsided, while galaxies with high concentration, mass, and density 
are not. We find that the strongest and most fundamental relationship between 
lopsidedness and the other structural parameters is with the surface mass den- 
sity. We also find, in agreement with previous studies, that lopsidedness tends to 
increase with radius. Both these results may be understood as a consequence of 
several factors. The outer regions of galaxies and low-density galaxies are more 
susceptible to tidal perturbations, and they also have longer dynamical times (so 
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lopsidedness will last longer). They are also more likely to be affected by any 
underlying asymmetries in the dark matter halo. 

Subject headings: galaxies: structure, galaxies: interactions, galaxies: general 



Introduction 



It has long been recognized that galaxies show large-scale asymmetries in their structure 
( iBaldwin et al.l Il980l ) . Lopsided galaxies have such asymmetries where one side of their 
disk is more massive and/or more extended tha n the opposite side. T his lopsidedness can 
be traced in the spatia l structure of the stars (IRix fc ZaritskyI Il995l ) and / or the HI gas 



( Richter fc Sancisilll994j ) and/or in the large-scale kinematics of this material (jSwaters et al. 
1999h . 



There are a variety of mechanisms or events that have been proposed to produce the 
observed lopsidedness. All of them involve a time-dependent non-equilibrium dynamical 
state, in most cases triggered through an external process. Such external processes are a 
natural consequence of the standard Lambda Cold Dark Matter cosmological framework. 
This implies that galaxies assemble hierarchically (a p rocess that is on-going). Examples 



that can lead to lopsidedness include a minor merger (jWalker et al.l Il996l : IZaritsky &: Rix 
19971 ). t he tidal interaction res ulting from a close encounter between roughly equal-mass 
galax ies (IKornreich et al.l 120021). and the asym metric accretion of intergalactic gas into the 



disk ( Bournaud et al. 



2005 



Keres et al.ll2005l ). Other mecha nisms involve the dark matter 



halo: stars and gas orbiting in a lopsided dark matter halo (IWeinbergill994l : I Jog] 119971 : |Jog 
1999h or a stellar/gas disk that is offset with res pect to the center of the dark matter halo 
( jLevine fc Sparke 1998 : Noordermeer et al. 2001). These also involve past tidal interactions 
and/or mergers that have perturbed the dark matter halo, but such perturbations may 
be quite long-lived. Finally, dynamical processes internal to the disk that lead to mildly 



lopsided distributions have also been investigated (IShu et al.l Il990l : ISyer fc Tremaind Il996 
Masset fc Taggei]ll997h . 



A variety of programs to study lopsidedness have been undertaken over the past decade. 
Most of these investigations have studied the lopsided dis tribution of the stellar compo- 
nent through analysis of optical and near-infrared images. IZaritsky fc Rixl (119971 ) studied 
a magnitude-limited sample of 60 field spiral galaxies. They measured lopsidedness as the 
radially averaged, azimuthal m = 1 Fourier amplitude Ai of the light (see Section [2^ below) 
and computed lopsidedness between 1.5 and 2.5 scale lengths in the galactic disks. The 
value of Al indicates the typical large-scale variation in mass density from side to opposite 
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side at the same distance from the galactic center. The mass density typically varies from 
between 1 ± times the average density at the same radius. They found tha t ~ 30% of field 



spiral galaxies exhibited significant lopsidedness {Ai > 0.2). iRudnick fc Rixl (119981 ) followed 
up this wo rk by studying lopside dness in 54 early-type galaxies and found that ~ 20% had 
Ai > 0.19. IConselice et al.l (120001 ) studied a sample of 113 z < 0.01 galaxies (elliptical, spiral, 
and irregular) but used a 180-degree rotational asymmetry measure Aisq. They found that 
asymmetry is strongly dependent on morphological type, with lower asymmetry in elliptical 
and lenti cular galaxies ari d high er asymmetry in late-type disk and irregular galaxies. More 
recently, iBournaud et al.l (120051 ) have measured the Fourier Ai parameter for 149 galaxies 
in the Ohio State University Bright Galaxy Survey. They confirmed that a large fraction 
of galaxies have significant lopsidedness in their stellar disks, with late-type galaxies being 
more lopsided. 

Lopsidedness in the light distribution can be produced by either a corresponding asym- 
metry in the underlying mass distribution in the stellar population or by large -scale variations 
i n the mass-to-light ratio (e.g., from star formation and dust obscuration). iRix fc Zaritsky 
(119951 ) investigated this issue with a sample of 18 face-on spiral galaxies imaged in the K' 
(2.2/i m) band where the effect of young stars or dust is minimized. They found that about 
a third of the sample s howed significan t lopsi dedness (similar to results from optical in- 
vestigations). Similarly, iRudnick fc Rixl (119981 ) found that lopsidedness in early-type disk 
galaxies is nearly identical when observed in the V, R, and I bands. They concluded that an 
asymmetric mass distribution then accounts for the majority of the asymmetry in the light 
distribution in these galaxies. 

Lopsidedness has also been studied in the distribution of HI gas. Since the HI can 
frequently be traced to significantly larger radii than the stars, these investigations are 
highly complementary to the optical image analysis. Due to the time-consuming nature 
of HI interferometri c mapping, only mo dest size samples have b een analyzed in this way 
(ISwaters et al.lll999l ). On the other hand, iRichter fc Sancisil (119941 ) have examined the global 
HI line profiles for roughly 1700 galaxies, and shown that at least 50% are significantly 
asymmetric (confirming that the large-scale HI distribution is frequently lopsided). HI maps 
also show th at apart from a lops ided distribution of the gas the HI rotation curves are often 
asymmetric (ISwaters et al.lll999l ). The connection between the phenona ena of structural and 
kinematic lopsidedness in galaxies is not yet clear (ISwaters et al.lll999l ). 



Despite these diverse investigations and the abundance of proposed models, the ori- 
gin of lopsidedness remains unsettled. For models involving tidal interactions or minor 
mergers, there is an expected link between lopsidedness and the loca l environment. The ev 



idence in this regard has been mixed (e.g., IWilcots fc PrescottI |2004J : iBournaud et al.l 12005 
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Andras et alJl2006l . I2OO7I : be Propris et all (120071 )1 



The investigations summarized above have all involved relatively small samples of galax- 
ies, making it difficult to assess the overall distribution of asymmetry or lopsidedness as a 
function of the basic parameters that characterize the structure of galaxies. This is the 
first of three papers in which we use the wealth of data available from the Sloan Digital 
Sky Survey (SDSS) to extend these studies of small samples (of-order one hundred galaxies) 
to large samples (tens of thousands). In this paper, we describe our sample selection and 
methodology. We also relate lopsidedness to the basic structural properties of the galaxies. 
In Paper II we will investigate the connection between lopsidedness and both star formation 
and black hole growth in galaxies. Finally, in Paper III we will examine the connection 
between lopsidedness and the local galaxy environment. 

In §21 we begin by presenting an initial low-redshift sample from the SDSS and de- 
scribe the observations and properties for its galaxies. Next, we explain our lopsidedness 
calculation. In ^|3l we address the major data quality issues that limit the reliability of the 
measurements for portions of the sample. On this basis, we apply cuts on the observational 
parameters to weed out the problematic cases for our subsequent analysis. Next, §1] describes 
the lopsidedness of galactic light distributions in different optical/near-IR bands, its corre- 
spondence with lopsided mass distributions, and its radial dependence. We then examine 
the relationship between lopsidedness and the basic structural properties of galaxies in §5] 
Finally, we summarize our findings in ^ 



2. Data & Basic Methodology 



2.1. Initial Sample 



The initial sample of ga laxies was taken from the Sloan Digital Sky Survey (lYork et al. 



2OOOI : IStoughton et al.ll2002l ). a large survey of photometric and spectro scopic data across tt 

sr. of the northern sky. The sample is derived fro m SDSS Data Rele ase 4 (lAdelman-McCarthy et al, 
20061 ). The survey's dedicated 2.5-n i telescope (Gunn et al.ll2006l ) at Apache Point Obser- 



vatory uses a unique CCD came ra (IGunn et al. 
r-, i-, and ^-band photometry (IFukugita et al 



19981 ) and drift-scanning to obtain u 



19961 : iHogg et al.l I2OOII : llvezic et al.l 12004 



Smith et al.l l2002l : iTucker et al.l 120061 ). The pixel scale is O'.'396/px. Fiber spectroscopy is 
obtained using 3" fibers and results in wavelength coverage between 3800—9200 A at a res- 
olution R = X/SX = 1850 - 2200. 

As we will show below, meaningful measurements of lopsidedness impose requirements 
on the angular size and signal-to- noise in the galaxy image. These criteria are not met 
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by the full SDSS Main galaxy sample (median redshift ~ 0.1 and magnitude r < 17.8; 
citealtswl+02). Accordingly, our initial galaxy sample was selected from the Main sample 
with a simple reds hift cut (z < 0-06). It contains 67107 galaxies. The SDSS photometric 
pipeline, PHOTO (ILupton et al.l 120011 ). provides Petrosian apparent magnitudes and half- 
and 90%-light radii {R^q and R^o) in each of the five bands, along with seeing conditions 
and other photometric a nd structural properties. We use the methodology described in 



Kauffmann et al.l (l2003al ) to use the spectral information to derive z-band mass-to-light 
ratios and hence stellar masses (M^). These stellar masses and the z-band half-light radii 
are then used to measure stellar surface mass densities (/i^, defin ed as the mean mass per 
unit area inside the half-light radius; see iKauffmann et al.ll2003bl ). 



2.2. Measuring Lopsidedness 



Lopsidedness in galactic light distribut ions has been compute d in recent years by two 
useful approaches. One approach, used by iRix fc ZaritskyI (Il995l ) and in later work, is to 
perform an azimuthal Fourier decomposition of galaxy light. The first mode, Ai, quantifies 
the large-scale overabundance of light in one side of a galaxy with a corresponding under- 
abundance in the opposite side. The first mode quantifies a galaxy-wide lopsidedness that 
can be quickly c ompu ted for large numbers of galaxies. In the second approach, detailed in 
Abraham et al.l (Il996l ). a different asymmetry index AigQ was devised by subtracting a 180° 
rotated image from the original galaxy image and summing the residual light compared to 
the total galaxy light. This lopsidedness measure is sensitive to both large- and small-scale 
variations in symmetry. To convert Aiqq into a measure of only large-scal e lopsidedness, 
a smoothing filter can be applied beforehand to the image (jConselicd l2003l ) . Our interest 
is to examine lopsidedness in tens of thousands of galaxies, and we have determined that 
the modal lopsidedness Ai algorithm is computationally less expensive than Aiso- We will 
therefore take the radially averaged first mode strength, Ai, as the measure of lopsidedness 
for the galaxy as a whole iii each band, as has been done in previous stu dies ( iRix fc Zaritsky 
I995I : IZaritskv fc Rb3 Il997l : iRudnick fc Rb3 Il998[ iRudnick et aPboOoh . The details of the 
calculation will fill out this section, but we first summarize the general method. 

The galaxy light is binned into radial and azimuthal bins, and then a finite Fourier 
transform is applied to the surface brightness fi in each radial bin. The transform for the 
kth radial bin is 

m=0 

The transform gives the mode magnitudes hmirk) and phases (pmiTk) for each mode m. 
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The magnitudes are divided by the zeroth mode (average surface brightness) and radially 
averaged to calculate for that galaxy and band, including the lopsidedness Ai. 

We compute this average over the radial range between R^q and Rqq. These are the 
practical limits: we can not go much inside R50 because of the seeing effects discussed below, 
and can not go much outside Rgo because the data become unacceptably noisy. As we will 
show below, the average radial variation in Ai beyond is modest, and the mean value for Ai 
does not depend significantly on the precise choices of the inner and outer radii. 

The transform also yields values of the higher-order modes. The second mode A2 rep- 
resents a combination of effects: ellipticity, inclination, two-arm spiral structure, and barred 
structure. Since we do not deproject the images to a face-on orientation, the amplitude of the 
second mode primarily measures the ellipticity of the image. Higher-order modes (m > 2) 
may include multi-arm spiral structure as well as Fourier "ringing" of features also present 
in lower-order modes. For example, a galaxy with a high ellipticity will have a strong second 
mode with weaker even-order modes also present. Similarly, a lopsided galaxy will have a 
strong first mode with weaker odd-order modes present as well. 

To begin our process, the galaxy center was determined. Precision centering is important 
because moving the center by more than 0.5 px can increase the lopsidedness significantly. 
The inner region of the galaxy was smoothed in a 3 px radius to reduce the effects of Poisson 
noise in determining the center. We assume that the center point of the galaxy is within 
3 px of the brightest pixel of the inner region of the galaxy. That brightest pixel is our 
initial estimate of the center point. To improve this estimate, the first moment of light was 
computed in a 3 x 3 px box centered on the brightest pixel. The first moment is the origin 
that minimizes the second moment and gives a central position in fractional units of a pixel 
that improves upon the brightest pixel estimate. 

Our centering method was used on each image without reference to the centers deter- 
mined for the same object in other bands. This means that each object was given a separate 
center point for each band. Nonetheless, the agreement between these independent center 
points is good. Of the three pairs of distances between centers in the three bands, the angular 
distance between the g- and i-band centers is the greatest. The la variation between these 



^ Conselice et al.l (|2000l ) have employed a different centering technique with their Aisa asymmetry measure. 
They found that the best center to use for the Aiso calculation was the center that minimizes Aigo- In our 
multi-mode calculation, we have found that the center that minimizes the strength of one mode does not 
minimize the strengths of other modes. Thus the minimization technique would require a different center 
for each mode. Since the modes are calculated simultaneously in the Fourier transform, minimization is not 
feasible. 



- 7- 



bands is 0'.'12too7 (0-32lo'i7 P^)- Some of the scatter may be attributable to the offset pixel 
locations between images of different bands, i.e., the coordinates of the center of the central 
pixel in one band will not be the located in the center of a pixel in another band. The rest 
of the scatter is due to variations in the central structure between the different bands. Our 
tests show that the mode strengths are significantly increased in a given image if the center 
point is moved more than 0.5 px from the best center. Thus the scatter of center points 
between the different bands should have no significant affect on the subsequent lopsidedness 
computation. 

The value of the sky background was calculated using the DAOPHOT package in the 
IDL Astronomy Library. For each galaxy and each band, the sky-subtracted image was 
partitioned into a polar grid of logarithmic radial bins between the Petrosian radii R^q 
and Rqo and centered at the determined galaxy center. We selected the minimum R^q and 
maximum Rgo determined in the g-, r-, and i-band images as the inner and outer radii of the 
grid, and the same grid size is used in the lopsidedness computation in all three bands. The 
radial bins were further partitioned into equal azimuthal bins. The number of radial bins 
was allowed to vary with the size of the galaxies by setting the innermost bin size as close to 
but no smaller than 1.0 pixel. Exactly 12 azimuthal bins were used in each grid so that the 
first six Fourier modes (m = 0, 1, 2, . . . , 5) could be determined. Each pixel that was split 
by the grid into more than one bin was divided into a 9 x 9 array of sub-pixels, and bilinear 
interpolation of neighboring pixels' surface brightness assigned each sub-pixel with a value 
of surface brightness. Higher-order interpolation, such as using bicubic splines, showed no 
significant improvement and thus was not worth the computational expense. The sub-pixels 
were then included in the appropriate polar bin. 

Galaxies whose observed light distributions are contaminated with light from other 
galaxies and foreground stars pose a problem for this Fourier decomposition. The contam- 
inating light results in a set of strong modes that describe the combination of galaxy and 
overlapping star or galaxy rather than the desired galaxy by itself. Foreground stars can 
be masked out with some ease once they arc located in images. We use the DAOPHOT 
routines of the IDL Astronomy Library to locate them and then estimate radii where their 
light drops to half the level of the la sky background noise. For intervening galaxies, we 
queried the SDSS PHOTO catalog for galaxy positions and Petrosian Rgg. Then a mask 
image was created that was the same size as the galaxy image, but each pixel accepted values 
of either (pixel is masked) or 1 (pixel is "good"). The Fourier decomposition calculation 
uses only the "good" pixels after the stars and galaxies are masked in the image. Any pixels 
that lay within a star, another galaxy, or other extraneous light source were ignored. 

After masking out the unwanted pixels, the surface brightness was computed in each 
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bin of the grid. An azimuthal Fourier transform (Eq. [T]) was performed in each radial bin 
to determine a magnitude bm{rk) and phase (pmi'^k) for each rn ode m. The transform was 



formulated as a general linear least squares problem (jPressll2002l ). Errors on the magnitudes 
and phases were deduced from the diagonal elements of the covariance matrix. The phase (pm 
of the mth mode is physically identical to phases (pm+'^T^K/m for integers K, so K was chosen 
for each radial bin to minimize the radial dependence of the phases. The magnitudes bm{rk) 
were then divided by the zeroth mode to give relative magnitudes amir^) = hm{rk) /hQijk)- 
These values arra('"fe) give a profile of lopsidedness (m = 1) and higher modes (m > 1), each 
normalized to the average brightness at radius r^. 

For simplicity, we prefer to use a single value of the Fourier modes in each band rather 
than a radial profile. We calculate the average of at each radius between R^q and i?9o, 
weighted by its error at each radius, to give a single, global measure (a^) for each galaxy 
in each band. The weighting naturally counts the dimmer, outer reaches of the galaxy less 
than the brighter, inner region and discounts any radial bin that is too contaminated with 
light from other sources. At radii where an extraneous light source has been masked out, 
the azimuthal Fourier transform can fail. Such radial bins are excluded from the weighted 
average. 

Since {am) is a positive definite quantity, random errors would preferentially overesti- 
mate the quantity. W e therefore correct the average in the manner used in the past (e.g.. 



Rudnick fc Rixlll998l ) and adopt 



Am = \l{amf-{6{am)f (2) 

as the strength of the mth mode, where 5 {am) is the random error in the weighted mean 
mode. We will discuss systematic and random errors in Ai in the next section (Q. 

As an example, we show the Fourier decomposition of the nearly face-on, barred spiral 
galaxy SDSS J125416.38-020204.4 in Fig. [H The first three panels show the g-, r-, and 
i— band images of this galaxy, a three-color image reconstructed from Fourier modes, and a 
three-color image of the m = 1 through 5 Fourier modes. The core of the galaxy is omitted 
in the three-color images as the decomposition is not attempted in that region. However, 
for this illustration, we have extended the Fourier decomposition to radii interior to R^q. 
The central bar is described by the strong even-order modes in all three bands at low radii 
2" < r < 5". The bar weakens at larger radii, and the even-order modes then decrease 
with increasing radii. At r > 8", the brightness of the galaxy drops until sky noise becomes 
significant. The lopsidedness in this galaxy is easily seen outside the bright center of the 
galaxy. The bar's brightness falls off more slowly with radius in the upper-right direction 
than in the lower-left direction. This lopsidedness is described in all three bands by the 
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significant first mode strength, which rises from near 0.1 to 0.4 between 2" and 8". The 
radially-averaged lopsidedness, as computed between R^q and Rqq in the i-band, is 0.29. 
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Fig. 1.— Images of the Fourier decomposition of SDSS J125416.38-020204.4 with A{ = 0.29. 
Upper row, left to right: SDSS g-, r-, and i-band images. Middle row, left: Combined gri 
image reconstructed from Fourier modes. Light from the ^f-band is colored blue, r-band 
green, and i-band red. Middle row, center and right: Combined gri images of the first and 
second modes alone. Lower row, left to right: Combined gri images of the third, fourth, and 
fifth modes alone. 
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Next, in Fig,[2l we show twelve example late- type galaxies along with their lopsidedness 
values in a progression from symmetric to lopsided. The three galaxies shown in the top 
row each exhibit a regular and symmetric appearance and have low values of lopsidedness 
Ai < 0.04. The next 5 galaxies have moderate Ai values of 0.08—0.13, and minor distortions 
are visible in each case. In cases of two-arm spiral galaxies, one or both arms are disturbed so 
that their shape and brightness are not identical at equal distances from the galactic center. 
In the case of SDSS J084845. 62-1-001729. 5 (third row, left), a second arm on the right has 
no clear counterpart on the left, and the lopsidedness of this galaxy is 0.11—0.12 depending 
on the band. The final four galaxies shown have clear asymmetry in brightness and shape 
and have Ai > 0.16 in each band. 

3. Final Sample Selection: The Effect of Systematic Errors 

Lopsidedness is ideally measured from an image of a bright, well resolved, and face-on 
galaxy without any overlapping background or foreground sources. In such an image, the 
Poisson noise of the detected light is negligible, and the size of the point-spread function is 
small compared to the characteristic size scale of the lopsidedness. Neither dust extinction 
nor inclination would significantly alter the light distribution. In real data, these conditions 
are not all met, and the resulting effects can systematically change the measurements of the 
Fourier modes. 

Below we address in turn how various sources of systematic error affect the strengths 
of Fourier modes. We then use this information to refine our sample selection and define a 
final sample of galaxies whose images are only negligibly affected by these systematic errors. 
These will form the basis of our scientific analysis. 

3.1. Ellipticity and Inclination 

We employ a circular polar grid to bin light from each galaxy, even if the galaxy appears 
elliptical and/or inclined on the image. The measured lopsidedness of elliptically projected 
galaxies can be underestimated by using the circular grid if the deprojected bright side of 
the galaxy coincides with the minor axis of the galaxy. This occurs because the increment 
and decrement of light along the minor axis due to lopsidedness is compared to the average 
light along a circular ring that also intersects the brighter, inner region at the major axis. 
Similarly, the lopsidedness will be overestimated if the bright and dim sides of the galaxy 
are aligned with the major axis. The resulting systematic error is diminished if the phase 



J101620.36 + 542847.1 
Al = 0.04. 0.04, 0.02 




J155547.64+423625.9 
Al = 0.08, 0.07, 0.06 




J08484S.62+001729.5 
Al = 0.12, 0.12, 0.11 



J124944.67-02S826.4 
Al = 0.20, 0.20, 0.22 
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J121935. 76+055048.3 
Al = 0.04, 0.03, 0.03 




J103705.40+124614.0 
Al = 0.08, 0.09, 0.08 




J085747.32-001159.9 
Al = 0.13, 0.13, 0.13 




|1-:4930.14+610654.2 
Al = 0.27, 0.27, 0.25 




J083104. 39+393721. 9 
Al = 0.03, 0.04, 0.04 



J225422. 30-101025.8 
Al = 0.11, 0.10, 0.10 




J232923. 56 + 141215. 4 
Al = 0.18, 0.16, 0.16 




J094706.75 + S41842.g 
Al = 0.26, 0.27, 0.26 




Fig. 2. — Twelve late-type SDSS galaxies shown with their g-, r-, and i-band lopsidedness. 
They are arranged by increasing lopsidedness toward the right in each row, least lopsided 
in the uppermost row and most lopsided in the lowermost row. Clear distortions in shape 
and/or brightness are visible in the galaxies with Ai > 0.08. 
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angle of the first Fourier mode wraps at least a quarter of the way around the galaxy. In 
this case, the first mode strength would be overestimated at some radii and underestimated 
at other radii, and the radially averaged Ai would have a smaller systematic error. 

Our sample has galaxies imaged in random orientations, and the over- and under- 
estimates of lopsidedness average to cancel each other except for highly inclined galaxies. 
We show this result in Fig. [31 where we have plotted the distribution of A\ (and also of A2) 
as a function of 6/a as measured in the z-band (similar results are seen in the g and r-bands). 

The relationship between b/a and A\ is clear. Face-on galaxies have a small second 
mode, typically 0.1-0.3, while highly inclined galaxies have a much higher second mode that 
often exceeds unity for b/a < 0.4. Though not shown, higher-order even modes show a 
similar trend with b/a as A2 but with weaker magnitude. 

Lopsidedness exhibits an increase at low b/a but is mainly independent of b/a for b/a > 
0.4. The rise of observed lopsidedness for increasingly inclined galaxies arises for a few 
reasons. First, dust lanes appear more optically thick when a galaxy is viewed nearly edge- 
on. A lopsided, edge-on distribution of dust may obscure the light in a symmetric galaxy 
to make the galaxy light appear lopsided. Second, the systematic error in using a circular 
grid for an elliptically projected galaxy becomes more severe at low b/a ratio. Nonetheless, 
A\ shows little dependence on b/a for b/a > 0.4, so our use of a circular polar grid and a 
finite-order transform on slightly elongated galaxies produces no significant systematic error 
when in reference to a large population of galaxies. The distributions of b/a as measured 
in the g-, r-, and z-bands are given in the lower panel of Fig. [31 A cut at 6/a = 0.4 in all 
three bands eliminates the 22% of the sample , and 52194 galaxies ar e retained. This cut on 



inclination is very similar to that adopted by iBournaud et al.l (l2005l ). 
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0.0 0.2 0.4 0.6 0.8 1.0 

b/a 



Fig. 3. — (Top two panels) Distributions of the Fourier modes Ai and A2 as functions of b/a 
in the i-band. The 5th, 25th, 50th, 75th, and 95th percentiles are shown. Lopsidedness is 
independent of inchnation for b/a > 0.4. On the other hand, A2, an indicator of elhpticity, 
inchnation, and two- arm spiral patterns, decreases significantly with b/a along the full range 
of b/a. (Lower right) The similar distributions of b/a as measured in the g- (green), r- 
(orange), and z-bands (red). A three-color cut at 6/a = 0.4 eliminates only a small portion 
of the sample. 
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Table 1. Example Galaxies in Figure [2] 



Row 


Cohimn Name A\ A\ 


1st 


Left 

Middle 

Right 


J101620.36+542847.1 0.04 0.04 0.02 
J121935.76+055048.3 0.04 0.03 0.03 
J083104.39+393721.9 0.03 0.04 0.04 


2nd 


Left 

Middle 

Right 


J155547.64+423625.9 0.08 0.07 0.06 
J103705.40+124614.0 0.08 0.09 0.08 
J225422.30-101025.8 0.11 0.10 0.10 


3rd 


Left 

Middle 

Right 


J084845.62+001729.5 0.12 0.12 0.11 
J085747.32-001159.9 0.13 0.13 0.13 
J232923. 56+141215.4 0.18 0.16 0.16 


4th 


Left 

Middle 

Right 


J124944.67-025826.4 0.20 0.20 0.22 
J144930. 14+610654.2 0.27 0.27 0.25 
J094706.75+541842.8 0.26 0.27 0.26 
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3.2. Spatial Resolution 

Lopsided galaxies systematically appear more symmetric when observed in conditions 
of poor seeing. Light is smeared from the brighter regions of the galaxy into surrounding, 
dimmer regions, reducing the contrast of the brighter and dimmer sides of the galaxy. The 
effect is more pronounced when the point-spread function is wide compared to the physical 
size of the galaxy. It is useful to define a seeing resolution 

^ 2i?50,& .ox 
PSF FWHMb ^ ' 

as the relative size of the galaxy compared to the FWHM of the PSF in band h. 

To demonstrate how poor seeing reduces lopsidedness, we have selected from our galaxy 
sample a subsample of 13500 isolated, well resolved galaxies with h/a > 0.5 and Si > 5.0 
and convolved the i-band images with a circular Gaussian PSF of 20 varying widths. We 
calculated A\ of each galaxy at each level of blurring. The subsample was binned by A\ of 
the original (not blurred) images. Fig. H] shows how the median A\ of the blurred images 
changes with seeing resolution in each bin of the original A\. Post-blurring A\ exhibits little 
change for highly resolved S > 5 galaxies. At moderate resolutions 2 < S* < 5, lopsidedness 
decreases with the seeing resolution. Barely resolved galaxies with S < 2 approach A\ ~ 0.05 
regardless of the lopsidedness at high resolution. 
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Fig. 4. — The i-band lopsidcdness of galaxies in images simulating various degrees of seeing 
resolution. A sample of isolated, wcll-rcsolved galaxies was divided into bins by Ai (measured 
before blurring) and the images were then blurred. For each bin, the relation between median 
Ai after blurring and the ratio of the half-hght and blurred PSF diameters (S) is shown 
above. Lopsided and symmetric galaxies can be distinguished for 5" > 3 but the range of Ai 
is diminished for -S" < 4. 
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The same effect is seen in our primary sample of galaxies without convolving the im- 
ages with a point-spread function. Fig. \5\ shows the distribution of A\ at different seeing 
resolutions along with the distributions of Sg, Sr, and Si for the whole sample (before the 
inclination cut is made). At low seeing resolution S < 3, the sample distribution of A\ shifts 
to lower lopsidedness (0.07 median), and few galaxies have A\ > 0.2. At high resolution 
S > 7, the distribution has a median that levels off around 0.12, and high lopsidedness 
A\ > 0.20 is found in 1 of every 4 galaxies. The 5th percentile increases with resolution to 
A\ = 0.05 at S = 7. 

Our simulations of poor seeing suggest that a distinction between symmetric and lop- 
sided galaxies can be identified by the Ai measure even if Ai is diminished by moderately 
poor seeing. The relative lopsidedness of well-resolved galaxies is preserved as the seeing is 
worsened down to S" ~ 3, i.e., the most symmetric galaxies with 3 < S* < 7 are likely to 
also be the most symmetric galaxies if they were observed instead at 5 > 7, and likewise 
for the more lopsided galaxies. However, the extremely lopsided galaxies may not be able 
to identified without higher resolution, perhaps S > 10. Trends at extreme lopsidedness 
Ai > 0.3 may not be reliable because only well resolved galaxies can lie in this lopsidedness 
range. At 5 = 4, seeing reduces the lopsidedness of 75% of galaxies by less than a factor 
of 1.5. We cut our sample at this seeing value and keep only galaxies better resolved than 
S = 4 in the g, r, and i bands. The cut alone reduces the sample size by ~ 37% to 42558 
galaxies. The more lopsided galaxies Ai > 0.2 are underrepresented due to this systematic 
effect. 
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Fig. 5. — Left: Distribution of i-band lopsidedness for galaxies observed at different seeing 
resolution. The 5th, 25th, 50th, 75th, and 95th percentiles are shown. Right: Distribution 
of seeing resolution as measured in the (7-band (green), r-band (orange), and i-band (red). 
Blurring from poor seeing systematically reduces the measured lopsidedness. A cut in see- 
ing resolution at 5* = 4 removes 37% of the sample while allowing a distinction between 
symmetric and asymmetric galaxies. 
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3.3. Flat-fielding error 

Poor flat-fielding in a galaxy image can increase the apparent lopsidedness of the galaxy. 
Poor correction of the CCD's spatial gradients in sensitivity can introduce a brighter sky 
level on one side of the galaxy and a dimmer sky level on the other side. Because a constant 
sky value is subtracted from the image in the lopsidedness computation, the gradient in sky 
brightness will increase the strength of the odd-order Fourier modes and overestimate the 
galaxy's true lopsidedness. We have calculated fiat-field errors in our images by calculating 
the sky level in regions in each of the four corners of the image and adopting the standard 
deviation of those sky levels as our fiat-fielding error. We find that fiat-fielding errors are 
typically negligibly small, ~ 1% of the sky value, and only become significant at large radii 
from a galactic center. In our calculation of Am, the contribution of the modes at large radii 
is weighted less, and so flat-fielding errors do not significantly affect the lopsidedness values 
that we have computed. We therefore do not impose a cut based on flat-fielding error for 
our SDSS-detected galaxies. 



3.4. Random Error 

Noisy data from low-surface-brightness galaxies can bias measurements of lopsidedness 
in a systematic way. The uncertainty in the correct light-weighted center will be higher 
in a dim galaxy than in a bright one, and a shifted center leads to an overestimate of the 
lopsidedness. 

We have undertaken two tests to determine the effect of noise on our measurements. 
First, we have compared the values of Ai measured in the r- and i-bands. Lopsidedness 
values have been dete rmined before to be mainly independent of wavelength in this range 



(iRudnick &: Rixlll998l ). and we will confirm this result for our sample as a whole in §4.1[ In 
the left panel of Fig. [6l we show the discrepancies between Ai in the two bands as a rough 
indicator of random error on the Ai measurement in general. Here the S/N ratio is computed 
for light within the i?5o-to-i?9o annulus in which Ai is computed. The interquartile (25th- 
75th percentile) range is small at high S/N (< 0.04 at S/N = 300) and gradually increases at 
low S/N as the distribution spreads out. At S/N = 30, the interquartile range has doubled 
to 0.09, and the median is skewed only slightly toward negative values of A^ — A\ ~ —0.01. 

We have applied another cut to our sample to remove dim, noisy galaxies with S/N < 30 
in any of the three gri bands. The right panel of Fig. E] shows the distribution of S/N in 
three bands. The cut removes a 5% of the sample. The sample retains 25155 galaxies after 
this cut and also the inclination and seeing cuts have been made. 
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Fig. 6. — (Left) Distributions of the difference in lopsidedness in the r and i bands as a 
function of i-band S/N. The 5th, 25th, 50th, 75th, and 95th percentiles are shown. The 
difference is typically small {\A[ — A\\ < 0.04) for most of the sample with S/Ni > 30. Since 
lopsidedness is systematically similar in these two bands, random errors in and A\ should 
be of the similar size as the difference | — A*]^ | . ( Right ) The distribution of S/ N in the g- 
(green), r- (orange), and i-band (red). A cut at S/N < 30 in all three bands eliminates a 
small fraction of the sample. 
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To further evaluate the effect of noise in calculating lopsidedness, we looked to the 
full SDSS DR4 dataset to retrieve a z < 0.06 sample of galaxies that have been observed 
3 or more times. We narrowed this sample down to a smaller sample where inclination, 
seeing resolution, and S/N met the same cuts as our main sample {b/a > 0.4, 5 > 4, 
and S/Ni^R^o < r < Rgo) > 30). We also required that the fractional RMS variations 
between different observations in seeing resolution and S/N were at most 10%. This latter 
requirement ensures that the repeated observations had similar observing conditions. The 
resulting sample contained 328 galaxies with repeated observations. We calculated the RMS 
differences in Ai in the repeated observations of these galaxies and adopted this measure as 
the error 6Ai in lopsidedness for these galaxies. 

Fig. [7| shows the distribution of 6Ai vs. Ai for the repeatedly observed galaxies. About 
a quarter of the galaxies show small errors in lopsidedness, 6Ai < 0.01, and median errors 
are typically ~ O.lAi. The largest errors are ~ O.SAi but affect only 5% of the sample. The 
third quartile rises as ~ 0.2y4i but can be as high as 0.02 at lopsidedness as low as 0.05. 
Typical errors can be expected to be the larger of 0.02 and O.lOAi. The majority of the 
sample has Ai < 0.20, and so 0.02 can be taken as the typical error for most of the sample, 
with larger errors for the most lopsided galaxies. This result is consistent with our estimate 
above based on comparing the r-band and i-band values for a much larger sample. 
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Fig. 7. — Errors in lopsidedness for galaxies with repeated observations. The 25th, 50th, 
and 75th percentiles are over-plotted. Typical lopsidedness errors are ~ 0.02 for the more 
symmetric galaxies and ~ 0.1 Ai for the more lopsided galaxies. 
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3.5. A Sample Suitable for Calculating Lopsidedness 

We have discussed the systematic errors that skew our calculation of lopsidedness. Here 
we summarize the proposed cuts to weed out cases where the systematic effects give unphys- 
ical values of Ai. After applying these cuts, we are left with a significant sample that is 
suitable for studying correlations between lopsidedness and other global properties of galax- 
ies. 



• Inclined galaxies. We eliminate all galaxies with b/a < 0.4 in any of the three bands 
to remove inaccurate Ai values computed using a circular aperture on an elliptically 
projected galaxy. This cut also removes cases of obscuration from optically thick dust 
lanes on edge-on galaxies. 

• Poorly resolved galaxies. We eliminate all galaxies with seeing resolution 5* < 4 in any 
of the three bands in order to remove galaxies whose lopsidedness is diminished due to 
poor seeing. 

• Dim galaxies. We eliminate all galaxies with S/N < 30 in the i?5o-to--R9o annulus 
(the region where Ai is calculated) in any of the three bands to remove cases where 
lopsidedness is augmented by Poisson noise and poor centering. 



Table [2] shows the census of galaxies retained after each cut is applied separately. The 
resolution cut alone removes the largest portion of the sample, and the noise cut removes the 
least amount. After all three cuts are applied to the initial sample of 67107 galaxies, 25155 
(37%) are retained. Unless otherwise specified, we have employed these cuts in our sample 
for all analysis presented below. These cuts depend on parameters linked to the observation 
of the galaxies but may also depend on physical properties. We next look at the structural 
properties of the sample and compare them to the larger DR4 sample from which it was 
drawn. 

The main structural parameters we will utilize in the analysis below are the stellar mass 
M^,, the effective stellar surface mass density fi^, (the mean stellar density interior to the 
z-band half light radius), and the con centration Ci (defined as the ratio of Rgo/R^Q in the 



i-band). See iKauffmann et al.l (l2003al ) for a detailed description. 



In Fig. [HI we show the distributions of these structural parameters M^, fi^, and Ci after 
we apply each cut in succession to the sample. The upper curve in each panel shows the 
distribution from the initial z < 0.06 DR4 sample. Then the inclination cut is applied (2nd 
distribution from the top), followed by the resolution cut (3rd distribution from the top) and 
the S/N cut (lower distribution). The inclination cut reduces the galaxy counts without any 
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strong correlation with structural properties. The resolution cut rejects more massive, high- 
mass-density, and concentrated galaxies, leaving the sample with a preference of late-type 
galaxies. Finally, the S/N cut causes little change in the relative proportions of massive, 
high-mass-density, or concentrated galaxies. 

The final sample contains 25155 galaxies spanning 3 orders of magnitude in stellar mass 
(10® — 10"*^^ Af^,), 3 orders of magnitude of stellar mass density (10^'^ — 10^'^ Mq kpc^^), and a 
wide range of z-band concentration (1.5 — 3.5), with a larger proportion of late- type galaxies 
than early-type galaxies. 



Table 2. Sample Cuts 



Cut 


Percentage 
Retained 


Galaxies 
Retained 


Initial Sample 


100% 


67107 


b/a > 0.4 


78% 


52194 


5 > 4 


63% 


42558 


S/N{R5o <r<Rgo)>30 


95% 


63434 


Final Sample after all 3 cuts 


37% 


25155 
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Fig. 8. — Distributions of the structural parameters Ci, M*, and /x* as the observational 
cuts were applied. From top to bottom, in each panel: full DR4 z < 0.06 sample before 
cuts were apphed; b/a > 0.4 cut is applied; S > 4 cut is applied; S/N > 30 cut is applied. 
The resulting sample has similar numbers of low- and high-mass galaxies, low- and high- 
stellar-density galaxies, but more a larger proportion of late-type galaxies than early-type 
galaxies. 
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4. Properties of Lopsidedness 

4.1. Light vs. Mass Distributions 

The Fourier modes describe the two-dimensional light distribution of galaxies. We wish 
to link these modes, specifically the lopsidedness, with a description of the two-dimensional 
mass distribution as seen along the same line of sight. If the contrast between the bright 
and dim halves of a galaxy is significantly different in the ^f-band than in the i-band, the 
lopsidedness measure may be indicating an asymmetry in mass-to-light ratios from asym- 
metrically distributed star formation and/or dust extinction. On the other hand, similar 
values in the two bands would suggest similar mass-to-light ratios and a corresponding lop- 
sidedness in surface mass density. To see which scenario is prevalent, we can look at the 
colors and magnitude differences between galaxies of different star formation histories and 
hence mass/light ratios. Then we can compare those differences to those of the bright and 
dim halves of the lopsided galaxies. 

We start with a set of pairs of stellar population models. We have taken 32000 simu- 
lations of stellar populat ions with identical mass and varying star formation histories from 



Kauffmann et al.l (j2003al ) and randomly paired them to look at their relative colors. The 
models were generated from a wide range of superimposed continuous and bursty star forma- 
tion histories with varying metallicity, such that bursty and continuously star-forming models 
each contribute about half of the models. We plot in the upper left panel of Fig. [9] the distri- 
butions of the difference in color A[g — i) and in magnitude Ai within the pairs. We find a 
tight relation between the difference in brightness and color within the pair {A{g—i) = 0.45Ai 
in the median). 

We next perform the same comparison using 113000 pairs of low-redshift {z < 0.06) 
observed galaxies. The pairs were selected from the SDSS DR4 galaxy sample and have 
been matched in redshift with Az < 0.001 and stellar mass with Alog^gM* < 0.01. The 
relation between relative color and magnitude is shown in the upper right panel of Fig. [91 
The brighter galaxy in the pair is the bluer galaxy, in agreement with the models. However 
the relation is somewhat shallower {A{g — i) = 0.29Az in the median). The difference in 
these relations may be due to a wider range of SFHs, especially the more extreme, bursty 
SFHs, in the models than in the observed galaxies. It may also reflect differences in the 
effects of dust extinction and reddening in the models vs. the data. 

We now compare these paired color-magnitude relations of the galaxies and the models 
to that determined from the bright and dim sides of the galaxies in our lopsided galaxy 
sample. Having determined the position angle that maximized the light asymmetry, we then 
measured the g, r, and i magnitudes for each galaxy half. Stars and intervening galaxies 
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were masked out and excluded, as before, but for each masked pixel here, the corresponding 
pixel at the same radius and 180° away is also excluded. Light is thus summed over the same 
area in each half. 

We show in the lower panel of Fig. [9] the difference in color A{g — i) between the con- 
trasting halves of galaxies in the DR4 sample as a function of their difference in i magnitude. 
The median color difference is zero and has little dependence on relative brightness. This is 
completely different from the behavior of the models or real galaxies as the mass/light ratio 
is va ried, and suggests that t he lopsidedness is not due to variations in mass/light ratio. In 



fact, iKauffmann et al.l (120071 ) have used SDSS galaxy spectra to show that the {g — i) fiber 
color is an excellent proxy for the stellar mass/light ratio. The fact that there is no system- 
atic difference in {g — i) color between the brighter and dimmer sides of the lopsided galaxies 
then implies that the lopsided light distribution is primarily tracing a lopsided stellar mass 
distribution. 

While there is no systematic offset in color with magnitude, we note that the spread in 
A{g~i) increases as Ai increases. This implies that there are large-scale spatial variations in 
the mass/light ratio in lopsided galaxies (presumably due to enhanced star formation and/or 
dust extinction). 
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Fig. 9. — Differences in g — i color and i magnitude for pairs of model stellar populations 
(upper left), for pairs of observed galaxies (upper right), and for bright and dim halves of 
observed galaxies (lower panel). The pairs of stellar populations and observed galaxies have 
random star formation histories and exhibit different mass-to-light ratios in the g and i 
bands. The contrasting halves of observed galaxies instead show no correlation between 
color and magnitude differences, implying similar average mass-to-light ratios in each half. 
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As another way of addressing the importance of variations in the the mass/hght ratio 
in causing lopsidedness, we have compared the distribution of Ai in three different bands 
in Fig. [TUJ The distribution measured from g-, r-, and i-band images are shown in green, 
orange, and red, respectively. 

The distributions of and A\ are nearly identical, but the Af distribution is skewed 
toward slightly higher values. We have already shown that the lopsidedness of the light 
distribution primarily traces the lopsidedness of the underlying mass distribution. Here we 
see the weaker, secondary effect. Newly formed stars are not uniformly distributed, and so 
the lopsidedness of the hght distribution includes a small contribution from the lopsidedness 
of the distribution of star formation. 




x 

ca 
CD 



CD 

E 




Fig. 10. — Distributions of the first and second azimuthal Fourier modes in the g— (green), 
r— (orange), and i— (red) bands between i?5o and Rqq. The r- and i-band distributions are 
similar. The ^f-band distribution is also similar but is skewed toward slightly higher values of 
lopsidedness. The lopsidedness of the light distributions is mainly tied to a lopsidedness in 
the mass distribution, but there is a small contribution from lopsidedness in the distribution 
of star formation. 
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4.2. Radial Dependence of Lopsidedness 



Previous studies h ave shown that asymmetry in galaxies shows a clear radial dependence. 
Rudnick fc Rixl (119981 ) used a sample of 54 face-on, ea rly-type disk galaxies and found that 
an increase in Ai lopsidedness with radius was normal. IConselice et al.l (120001) studied radial 
profiles of both elliptical and disk galaxies in a sample of 113 from the iFrei et al.l (119961 ) 
sample. They found that asymmetry typically increases at larger radii r > R^o in disk 
galaxies but peaks at smaller radii in elliptical and lenticular galaxies. One would expect 
our larger sample of galaxies to confirm the dependence of lopsidedness on both radius and 
Hubble type. In what follows, we use concentration as a measure of Hubble type. 
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Fig. 11. — Radial profiles of lopsidedness shown as distributions of function of 

radius. The first three panels show distributions for symmetric, average, and lopsided late- 
type galaxies (Cj < 2.6), and the lower right panel shows the distribution for all early-type 
galaxies (Cj > 2.6). The 5th, 25th, 50th, 75th, and 95th percentiles are shown. Lopsidedness 
increases with radius in symmetric and lopsided, late-type and early-type galaxies alike. 
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In Figim we show radial profiles as the distribution of the first Fourier mode (lowercase) 
a\{r) at radii extending outward to Rgo. The first three panels show profiles for late-type 
galaxies (Cj < 2.6) partitioned into bins of global lopsidedness (capitalized A\), and the last 
panel shows early-type (Cj > 2.6) galaxy profile. R^q corresponds to roughly 0.38i?9o to 
0.60i?9o for these late-type galaxies. In all these profiles, lopsidedness increases steadily at 
radii larger than ~ O.S-Rgo. Moderately and highly lopsided {A\ > 0.08) late-type galaxies 
also show a gradual increase in lopsidedness with radius at smaller radii. The early-type 
and more symmetric late-type galaxies show a minor decline in lopsidedness at small radii 
and an increase at larger radii. The decline at small radii may be an effect originating from 
minor centering errors near the galactic center, causing a small overestimation of a\{r) at 
the smallest radii. This effect diminishes with radius. At large radii, the increase in a\{r) 
is tied to a real increase in light asymmetry, though there could also be some unphysical 
enhancement due to the low S/N. 



5. Structural Properties of Lopsided Galaxies 



In §3.51 we showed that our working sample was drawn from a full sample of 2; < 0.06 
galaxies and was selected based on cuts on several observational parameters. The proportion 
of high-mass-density and highly concentrated galaxies was reduced, but a significant number 
of these galaxies was also retained. Our sample thus allows us to study lopsidedness of 
galaxies over a wide range of the basic structural properties of the galaxies, namely their 
concentration, stellar mass, and stellar mass density. 

The concentration is a rough proxy for Hubble type, with higher values corresponding 
to earlier type s. The corre spondence between d a nd Hubble Type has been considered by 
Strateva et al.l (1200 ll ) and IShimasaku et al.l (I2OOII ). The correspondence is not tight, but 
the value Cj = 2.6 is th e rough dividing line between early- and late-type galaxies (see also 
Kauffmann et al.ll2003al . The top- left panel of Fig. [T^ shows the distribution of the global A\ 
as a function of concentration. At Cj = 2.6, A\ has a moderate value near 0.10. For Cj > 3.1 
(typical of elliptical galaxies), the lopsidedness is typically too small to reliably measure. 
Below Cj = 2.6 the lopsidedness rises systematically with decreasing concentration (later 
Hubble types). Lopsidedness is thus commonplace in the late-type field galaxies and the 
galaxies in low-density environments that dominate o ur sample, in agreeme nt with previous 



studies ( iMatthews et al.lll998l : IConselice et al.ll2000l: iBournaud et a. 



result might not extend to other environments (lAngiras et al. 
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a 



20071 ). 



2OO5I ). However, this 



The top -right panel of Fig. [T2] sh ows the distribution of lopsidedness as a function of 
stellar mass. iKauffmann et al.l (l2003bl ) showed that the local galaxy population is bimodal, 
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with the transition from late type star forming galaxies to early type old galaxies occurring 
at a value of ~ IO^'^'^Mq. We see that lopsidedness is only significant for the low mass 
population, and increases systematically with decreasing mass. 

Finally, we show that there is a similar but even strong er trend between l o psided ness and 
stellar mass density, as shown in the lower panel of Fig. [T2j iKauffmann et al.l (l2003bl ) showed 
that the transition between late type star forming galaxies and early type old galaxies occurs 
at a value /i* ~ 10^'^ Mq kpc~^. We find that lopsidedness is only commonplace among the 
low density population. The lopsidedness increases very strongly with decreasing density 
in this regime. Indeed at the lowest densities (yU* < lO"^'^ Mq kpc~^) typical galaxies are 
significantly lopsided {A\ > 0.20), and few galaxies in this density range are less lopsided 
than 0.10. 
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Fig. 12. — Distributions of i-band lopsidedness as functions of the structural parameters 
M*, and /x*. The 5th, 25th, 50th, 75th, and 95th percentiles are shown in each panel. 
Lopsided light distributions are more common in low-mass, low-mass-dcnsity, and low- 
concentration galaxies. Massive, dense, and concentrated galaxies tend to be considerably 
more symmetric. 
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The above st ructural parameters of SP SS galaxies (C, M*, /i*) are highly correlated with 



one another (e.g.. iKauffmann et al.ll2003af ). From the above plots, it is not immediately clear 
whether there is a separate underlying physical correlation between lopsidedness and each 
of the structural parameters, or whether some of the apparent correlations are induced by 
mutual correlations between the parameters. To explore this, we now turn to Fig. [13] and 
compare the three structural properties, two at a time, with lopsidedness. First, we have 
plotted /i* and Ci against in two-dimensional bins and color-coded each bin holding 
at least 5 galaxies by its median A[. We see immediately that while lopsidedness does 
indeed simultaneously correlate with all these structural parameters, the correlation with yU* 
is the strongest and most fundamental. In the plots of /i* vs. and Ci we see that the 
lopsidedness at a given value of is essentially independent of either M^, or Q. Conversely, 
for given values of M^, or Ci, there is a systematic increase in lopsidedness with decreasing 
fi^. The plot of Ci vs. implies that both parameters are correlated with lopsidedness 
(neither is clearly the more fundamental). 

We have used linear partial correlation analysis to see which pairs of these four structural 
properties show more fundamental correlations. We calculated the linear (Pearson) corre- 
lation coefficients for each of the six combinations of the four parameters. We use logM^, 
log/i*, \ogA\, and Ci for this analysis because the relationships are more linear if logarithms 
are used for three of the parameters. For each pair, we also removed the dependence of the 
remaining two parameters and calculated the partial correlation coefficient. The coefficients 
are listed in Table [31 Lopsidedness correlates moderately with all three of the other param- 
eters but most strongly with stellar mass density (corr. coeff. = —0.56), then concentration 
(—0.47) and mass (—0.46). Once the correlations with stellar mass and concentration are 
removed, the correlation between lopsidedness and mass density is reduced in magnitude to 
—0.20. The partial correlations between lopsidedness and the other two structural param- 
eters are weaker (0.00 with Ci and —0.12 with log M^,). In agreement with Fig. [T31 we see 
that the most fundamentral correlation between lopsidedness and a structural parameter is 
with the surface mass density. 
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Fig. 13. — Stellar mass-mass density and concentration-mass relationships colored by median 
i— band lopsidedness. Galaxies that are massive, dense, and concentrated tend to have sym- 
metric light distributions, while low-mass, low-mass- density, and low- concentrated galaxies 
tend to be lopsided. 
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Table 3. Correlation Coefficients of Structural Parameters 



Dependence (Partial) 



Par. 1 


Par. 2 




Removed 


Corr. Coeff. 


logA\ 


log/x* 






-0.56 


log 4 


log/U* 


log 




-0.20 


log A\ 


logM* 






-0.46 


log 4 


logM^ 


log 


in* k Ci 


0.00 


log 4 


Ci 






-0.47 


log 4 


Ci 


log 


; /X* &; log M* 


-0.12 


logyU* 


logM* 






0.87 


log/i* 


logM* 


log 




0.82 


log/i* 








0.70 


log/i* 




log 


lA[ k logM* 


0.59 


logM* 








0.49 


logM* 




log 


; A\ & log /X* 


-0.32 
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6. Conclusions &z Discussion 

We have measured large-scale galactic asymmetry for a large sample of low-redsliift {z < 
0.06) galaxies drawn from the Sloan Digital Sky Survey. Our use of lopsidedness, a radially 
averaged m = 1 azimuthal Fourier mode, has proven useful for a large fraction of the sample. 
Images of a minority of galaxies in the sample have poor observational properties that cause 
significant systematic errors in the lopsidedness calculation, and these galaxies were removed 
from the sample via cuts on angular size, signal-to- noise, and ellipticity/inclination. Those 
cuts removed a higher fraction of the high-mass, high-mass-density, and high-concentration 
galaxies than those with low values of these structural properties. Nonetheless, the resulting 
sample is well represented by the galaxies of the same range of structural properties as the 
original sample. 

We find that there are no systematic differences between the (g — i) colors of the brighter 
and fainter sides of lopsided galaxies. Th is implies that there is no systematic difference in the 



mass/light ratio (IKauffmann et al.ll2007l ). and hence that the lopsided light distributions are 
primarily caused by lopsided distributions in the stellar mass. We have verified this through 
analysis of the relationship between color and mass/light ratio for both model galaxy spectral 
energy distributions and SDSS galaxy data. However, for our sample the lopsidedness in the 
(jf-band tends to be slightly greater than in the r- and z-bands. Thus, some of the lopsidedness 
in the light does arise from the effects of star-formation and/or dust extinction (which will 
more strongly affect the (yf— band light). 

Lopsidedness is a structural property that depends strongly on other structural prop- 
erties. Galaxies with progressively lower concentration, stellar mass, or stellar mass density 
tend to have progressively higher lopsidedness. We show that the strongest and most fun- 
damental correlation is between lopsidedness and stellar mass density. We also find that 
lopsidedness increases systematically with increasing radius, particularly for late-type galax- 
ies. 

Lopsidedness can be induced through tidal stress associated with interactions with 



a companion galaxy or through accretion or minor mergers (e.g. IZaritsky &: Rixl Il997 



Bournaud et al.l l2005l ). Galaxies with low density will be most affected by tidal stress, 
and the effects of a tidal perturbation will last longer in such systems due to the longer 
dynamical times. The same arguments pertain to the outer parts of galaxies. Thus, the two 
above results make good physical sense. Alternatively, if the dark matter halo is lopsided, its 
effects on the structure of the stellar disk will be more pronounced in the outer region and 
in galaxies with low mass and low density (where dark matter is more dynamically impor- 
tant). The relatively large values of lopsidedness we measure to be commonplace {Ai > 0.1) 
appear to be too large to be generated by internally generated dynamical processes (e.g.. 
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Masset fc Taggei]ll997h . 



Our overall goal in this investigation has been to use lopsidedness as a way of quantifying 
the signature of moderate or weak global dynamical perturbations. The next step will be 
to determine the connections between such perturbations and both the on- going/recent star 
formation and the growth of supermassive black holes in galaxies. These connections can help 
constrain the processes and conditions that guide the formation and evolution of the galaxies. 
In future papers we will address these questions using the present sample of galaxies. 
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